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ABSTRACT 

• i— i , We constrain flat cosmological models with a joint likelihood analysis of a new compila- 

■ tion of data from the cosmic microwave background (CMB) and from the 2dF Galaxy 
5h I Redshift Survey (2dFGRS). Fitting the CMB alone yields a known degeneracy be- 
tween the Hubble constant h and the matter density f2 m , which arises mainly from 
preserving the location of the peaks in the angular power spectrum. This 'horizon- 
angle degeneracy' is considered in some detail and shown to follow a simple relation 
n m h 3A = constant. Adding the 2dFGRS power spectrum constrains Q m h and breaks 
the degeneracy. If tensor anisotropies are assumed to be negligible, we obtain values 
for the Hubble constant h = 0.665 ±0.047, the matter density fl m = 0.313 ±0.055, and 
the physical CDM and baryon densities fl c h 2 = 0.115 ± 0.009, Q. b h 2 = 0.022 ± 0.002 
(standard rms errors). Including a possible tensor component causes very little change 
to these figures; we set an upper limit to the tensor-to-scalar ratio of r < 0.7 at 95% 
confidence. We then show how these data can be used to constrain the equation of state 
of the vacuum, and find w < —0.52 at 95% confidence. The preferred cosmological 
model is thus very well specified, and we discuss the precision with which future CMB 
data can be predicted, given the model assumptions. The 2dFGRS power-spectrum 
data and covariancc matrix, and the CMB data compilation used here, are available 



from tittp://www. roe. ac.uk/~wjp/ 
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1 INTRODUCTION 

The 2dF Galaxy Redshift Survey (2dFGRS; see e.g. Col- 
less et al. 2001) has mapped the local Universe in detail. If 
the galaxy distribution has Gaussian statistics and the bias 
factor is independent of scale, then the galaxy power spec- 
trum should contain all of the available information about 
the seed perturbations of cosmological structure: it is statis- 
tically complete in the linear regime. The power spectrum 
of the data as of early 2001 was presented in Percival et al. 
(2001), and was shown to be consistent with recent cosmic 
microwave background (CMB) and nucleosynthesis results. 

In Efstathiou et al. (2002) we combined the 2dFGRS 
power spectrum with recent CMB datasets in order to con- 
strain the cosmological model (see also subsequent work by 
Lewis & Bridle 2002). Considering a wide range of possible 
assumptions, we were able to show that the universe must 
be nearly flat, requiring a non-zero cosmological constant A. 
The flatness constraint was quite precise (|1 — fitot| < 0.05 
at 95% confidence); since inflation models usually predict 
near-exact flatness (|1 — fitot| < 0.001; e.g. Section 8.3 of 
Kolb & Turner 1990) , there is strong empirical and theoret- 
ical motivation for considering only the class of exactly flat 
cosmological models. The question of which flat universes 
match the data is thus an important one to be able to an- 
swer. Removing spatial curvature as a degree of freedom also 
has the practical advantage that the space of cosmological 
models can be explored in much greater detail. Therefore, 
throughout this work we assume a universe with baryons, 
CDM and vacuum energy summing to fi t ot = 1 (cf. Peebles 
1984; Efstathiou, Sutherland & Maddox 1990). 

In this work we also assume that the initial fluctua- 
tions were adiabatic, Gaussian and well described by power 
law spectra. We consider models with and without a tensor 
component, which is allowed to have slope and amplitude 
independent of the scalar component. Recent Sudbury Neu- 
trino Observatory (SNO) measurements (Ahmad et al. 2002) 
are most naturally interpreted in terms of three neutrinos 
of cosmologically negligible mass (< 0.05 eV, as opposed 
to current cosmological limits of order 2 eV - see Elgaroy 
et al. 2002). We therefore assume zero neutrino mass in this 
analysis. In most cases, we assume the vacuum energy to 
be a 'pure' cosmological constant with equation of state 
w = p/pc 2 = — 1, except in Section ^ where we explore 
io > -1. 

In Section ^ we use a compilation of recent CMB ob- 
servations (including data from VSA (Scott et al. 2002) and 
CBI (Pearson et al. 2002) experiments) to determine the 
maximum-likelihood amplitude of the CMB angular power 
spectrum on a convenient grid, taking into account calibra- 
tion and beam uncertainties where appropriate. This com- 
pression of the data is designed to speed the analysis pre- 
sented here, but it should be of interest to the community 
in general. 

In Section |^ we fit to both the CMB data alone, and 
CMB + 2dFGRS. Fits to CMB data alone reveal two well- 
known primary degeneracies. For models including a possi- 
ble tensor component, there is the tensor degeneracy (Efs- 
tathiou 2002) between increasing tensors, blue tilt, increased 
baryon density and lower CDM density. For both scalar-only 
and with-tensor models, there is a degeneracy related to the 
geometrical degeneracy present when non-flat models are 



Table 1. Best-fit relative power calibration corrections for the 
experiments considered are compared to expected rms errors. In 
addition, we recover a best fit beam error for BOOMERaNG of 
+0.4%, measured relative to the first data point in the set, and 
+0.07% for Maxima. 



power calibration error 
experiment best-fit (%) rms (%) 



BOOMERaNG 


-13.5 


20 


Maxima 


+1.6 


8 


DASI 


+0.9 


8 


VSA 


-0.3 


7 


CBI 


+0.7 


10 



considered, arising from models with similar observed CMB 
peak locations (cf. Efstathiou & Bond 1999). In Section ^ 
we discuss this degeneracy further and explain how it may 
be easily understood via the horizon angle, and described 
by the simple relation Q. m h iA — constant. 

Section [| considers a possible extension of our standard 
cosmological model allowing the equation of state param- 
eter w of the vacuum energy component to vary. By com- 
bining the CMB data, the 2dFGRS data, and an external 
constraint on the Hubble constant h, we are able to con- 
strain w. Finally, in Section ^, we discuss the range of CMB 
angular power spectral values allowed by the present CMB 
and 2dFGRS data within the standard class of flat models. 



2 THE CMB DATA 

Recent key additions to the field of CMB observations come 
from the VSA (Scott et al. 2002), which boasts a smaller 
calibration error than previous experiments, and the CBI 
(Pearson et al. 2002, Mason et al. 2002), which has extended 
observations to smaller angles (larger £'s). These data sets 
add to results from BOOMERaNG (Netterfield et al. 2002), 
Maxima (Lee et al. 2001) and DASI (Halverson et al. 2002), 
amongst others. Rather than compare models to each of 
these data sets individually, it is expedient to combine the 
data prior to analysis. This combination often has the ad- 
vantage of allowing a consistency check between the indi- 
vidual data sets (e.g. Wang et al. 2002). However, care must 
be taken to ensure that additional biases are not introduced 
into the compressed data set, and that no important infor- 
mation is lost. 

In the following we consider COBE, BOOMERaNG, 
Maxima, DASI, VSA and CBI data sets. The 
BOOMERaNG data of Netterfield et al. (2002) and the 
Maxima data of Lee et al. (2002) were used assuming the 
data points were independent, and have window functions 
well described by top-hats. The £ < 2000 CBI mosaic field 
data were used assuming that the only significant corre- 
lations arise between neighbouring points which are anti- 
correlated at the 16% level as discussed in Pearson et al. 
(2002). Window functions for these data were assumed to 
be Gaussian with small negative side lobes extending into 
neighbouring bins approximately matched to figure 11 of 
Pearson et al. (2002). We also consider the VSA data of 
Scott et al. (2002), the DASI data of Halverson et al. (2002), 
and the COBE data compilation of Tegmark et al. (1996), 
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Figure 1. Top panel: the compilation of recent CMB data used in our analysis (see text for details). The solid line shows the result of a 
maximum-likelihood fit to these data allowing for calibration and beam uncertainty errors in addition to intrinsic errors. Each observed 
data set has been shifted by the appropriate best-fit calibration and beam correction. Bottom panel: the solid line again shows our 
maximum-likelihood fit to the CMB power spectrum now showing the nodes (the points at which the amplitude of the power spectrum 
was estimated) with approximate errors calculated from the diagonal elements of the covariance matrix (solid squares). These data are 
compared with the compilation of Wang et al. (2002) (stars) and the result of convolving our best fit power with the window function 
of Wang et al. (crosses). In order to show the important features in the CMB angular power spectrum plots we present in this paper we 
have chosen to scale the x-axis by (log£) 2,5 . 



for which the window functions and covariance matrices 
are known, where appropriate. The calibration uncertain- 
ties used are presented in Table [jj and the data sets are 
shown in Fig. nj. In total, there are 6 datasets, containing 68 
power measurements. 

In order to combine these datasets, we have fitted a 
model for the true underlying CMB power spectrum, con- 
sisting of power values at a number of I values (or nodes). 
Between these nodes we interpolate the model power spec- 
trum using a smooth Spline3 algorithm (Press et al. 1992). 
The assumption of smoothness is justified because we aim 
to compare CMB data with CDM models calculated us- 
ing CMBFAST (Seljak & Zaldarriaga 1996). Internally, this 
code evaluates the CMB power spectrum only at a partic- 
ular set of I values, which are subsequently Spline3 inter- 
polated to cover all multipoles. It is therefore convenient to 
use as our parameters the CMB power values at the same 



nodes used by CMBFAST in the key regime 150 < I < 1000. 
By using the same smoothing algorithm and nodes for our 
estimate of the true power spectrum, we ensure that no 
additional assumptions are made in the data compilation 
compared with the models to be tested. For I < 150 and 
£ > 1000 the data points are rather sparsely distributed, 
and we only selected a few £ values at which to estimate the 
power. The best-fit amplitude of the power spectrum at an 
extra node at £ — 2000 was determined in our fit to the ob- 
served CMB data, in order that the shape of the interpolated 
curve around I = 1500 had the correct form. This was sub- 
sequently removed from the analysis, and models and data 
were only compared for £ < 1500. In addition to requiring no 
interpolation in CMBFAST, this method of compression has 
a key advantage for our analysis. Normally, CMB data are 
expressed as bandpowers, in which one specifies the result 
of convolving the CMB power spectrum with some kernel. 
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Table 2. Recovered best-fit power spectrum values with rms 
values given the 6 data sets analysed. 





ST 2 / [iK 2 


rms error / 


2 


314 


443 


1 


803 


226 


8 


770 


156 


15 


852 


174 


50 


1186 


1414 


90 


2796 


673 


150 


3784 


546 


200 


5150 


627 


250 


5306 


590 


300 


3407 


364 


350 


2339 


265 


400 


1627 


205 


450 


1873 


202 


500 


2214 


240 


550 


2479 


249 


600 


2061 


245 


650 


1849 


244 


700 


2023 


274 


750 


1614 


295 


800 


2089 


373 


850 


2654 


475 


900 


2305 


515 


950 


1178 


480 


1000 


1048 


320 


1200 


1008 


214 


1500 


530 


178 



This remains true of some previous CMB data compilations 
(e.g. Wang et al. 2002). In contrast, we estimate the true 
power spectrum at a given I directly, so that no convolution 
step is required. This means that parameter space can be 
explored more quickly and efficiently. 

Given a set of nodal values, we form an interpolated 
model power spectrum, convolve with the window function 
of each observed data point and maximized the Likelihood 
with respect to the nodal values (assuming Gaussianity - 
see Bond, Jaffe & Knox 2000 for a discussion of the possible 
effect of this approximation). Calibration errors and beam 
uncertainties were treated as additional independent Gaus- 
sian parameters, and were combined into the final likelihood, 
as well as being used to correct the data. The resulting best 
fit calibration and beam errors are compared to the expected 
rms values in Table [[]. 

In agreement with Wang et al. (2002), we find a neg- 
ative best fit BOOMERaNG calibration correction (13% 
in power), caused by matching data sets in the regime 
300 < I < 500. Applying this correction (included in the 
data points in Fig. ^) slightly decreases the amplitude of 
the first peak. Nevertheless, our combined power values are 
systematically higher than in the compilation of Wang et 
al. (see the lower panel of Fig. 1). This derives partly from 
the inclusion of extra data, but also results from a bias in 
the analysis method of Wang et al. They use the observed 
power values to estimate the error in the data, rather than 
the true power at that multipole (which we estimate from 
our model). A low observational point is thus given a spuri- 
ously low error, and this is capable of biasing the averaged 
data to low values. 

The final best fit power spectrum amplitudes given 



the 6 data sets analysed are presented in Table |2j with 
the corresponding ^-values of the nodes and rms errors. 
Formally this fit gave Xmin — 31.9, given 34 degrees of 
freedom (there are 68 data points, and we estimate 27 
power spectrum values, 5 calibration and 2 beam correc- 
tions) . This result demonstrates that the different data sets 
are broadly consistent after allowing for calibration and 
beam uncertainty. The Hessian matrix of the likelihood 
provides an estimate of the inverse covariance matrix for 
the power spectrum estimates. This was calculated numer- 
ically and_Js_avajJabl£ J _togeUier_with the averaged data, 
from ittp: //www .roe . ac.uk/~wjp/ As emphasised previ- 
ously, these are estimates of the true power at the £ values 
given and therefore do not require window functions. In the 
following Section we use these CMB results to constrain flat 
cosmological models. 



3 COSMOLOGICAL MODELS 
3.1 Parameter space 

In the following we parametrise flat cosmological models 
with seven parameters (plus two amplitudes): these are the 
physical baryon densityrl Q,bh 2 , the physical CDM density 
Q c h 2 , the Hubble constant h, the optical depth to the last 
scattering surface r, the scalar spectral index n s , the ten- 
sor spectral index nt and the tensor-to-scalar ratio r. The 
tensor-to-scalar ratio r is defined as in Efstathiou et al. 
(2002): the scalar and tensor GYs are normalized so that 



1UUU 

h 5^ (2f + 1)( ^ s = (4 x 10 ~ 5)5 

i- f> + 1)CJ = (2xl0- 5 ) 2 



(1) 
(2) 



Ce is then given by Ce = Q 2 (Cf + rCj), where Q 2 is the 
normalization constant. We marginalize over both this and 
the amplitude of the 2dFGRS power spectra in order to 
avoid complications caused by galaxy biasing and redshift 
space distortions (Lahav et al. 2002). 

CMB angular power spectra have been calculated us- 
ing CMBFAST (Seljak & Zaldarriaga 1996) for a grid of 
~ 2 x 10 s models. For ease of use, a uniform grid was adopted 
with a varying resolution in each of the parameters (details 
of this grid are presented in Table []) . Likelihoods were cal- 
culated by fitting these models to the reduced CMB data set 
presented in Section |[ Similarly, large scale structure (LSS) 
power spectra were calculated for the relevant models us- 
ing the fitting formula of Eisenstein & Hu (1998), and were 
convolved with the window function of the 2dFGRS sample, 
before being compared to the 2dFGRS data as in Percival 
et al. (2001). 

In order to constrain parameters, we wish to determine 
the probability of each model in our grid given the available 



* As usual, f2;,,n c are the densities of baryons & CDM in units 
of the critical density, and h is the Hubble constant in units of 
lOOkms - 1 Mpc - . 'Derived' parameters include the matter den- 
sity n m = Q c + fi;,, and Q\ = 1 — fi ra . 
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Figure 2. Two parameter likelihood surfaces for scalar only models. Contours correspond to changes in the likelihood from the maximum 
of 2A In C = 2.3, 6.0, 9.2. Dashed contours are calculated by only fitting to the CMB data, solid contours by jointly fitting the CMB and 
2dFGRS data. Dotted lines show the extent of the grid used to calculate the likelihoods. 



Table 3. The distribution of parameters (defined in the text) in 
the ~ 2 X 10 s flat cosmological models considered in this paper. 
The grid used was linear in each parameter between the limits 
given in order to simplify the marginalization assuming a uniform 
prior on each. 



parameter 


min 


max 


grid size 




0.01 


0.04 


25 


tt c h 2 


0.05 


0.22 


25 


h 


0.40 


1.00 


25 


T 


0.00 


0.10 


2 


n s 


0.80 


1.30 


25 


nt 


-0.20 


0.30 


10 


r 


0.00 


1.00 


25 



CMB and 2dFGRS data. However, we can only easily calcu- 
late the probability of the data given each model. In order 
to convert between these probabilities using Bayes theorem, 
we need to adopt a prior probability for each model or pa- 
rameter. In this work, we adopt a uniform prior for the pa- 
rameters discussed above between the limits in Table |3| i.e. 
we assume that the prior probability of each model in the 
grid is the same. Assuming a uniform prior for physically 
motivated parameters is common in the field, although not 
often explicitly mentioned. Note that the constraints placed 
by the current data are tight compared with the prior, and 
that the biases induced by this choice are therefore relatively 
small. 

The likelihood distribution for a single parameter, or 



for two parameters can be calculated by marginalizing the 
estimated probability of the model given the data over all 
other parameters. Because of the grid adopted in this work, 
we can do this marginalization by simply averaging the C 
values calculated at each point in the grid. 

In Fig. |^ we present two-parameter likelihood contour 
plots (marginalized over the remaining parameters) for the 
subset of scalar-only models i.e. r fixed at 0. For these scalar- 
only models, we choose to plot r only against f2(,ft 2 as r is 
poorly constrained by the CMB data, and has no degenera- 
cies with the other parameters. In Fig. ^ we present two- 
parameter likelihood contour plots (marginalized over the 
remaining parameters) for models allowing a tensor compo- 
nent. The spectral index of the tensor contribution is poorly 
constrained by the CMB data so, as for r, we only show one 
plot with this parameter. 

Figs ^ & ^ reveal two key directions in parameter space 
that the CMB data have difficulty constraining. When a ten- 
sor component is included, we have the tensor degeneracy - a 
trade-off between increasing tensors, increasing n s , increas- 
ing Qbh 2 and decreasing Q c h 2 (for more detail see Efstathiou 
2002). In addition, in both the scalar-only and with-tensor 
cases, there is a degeneracy between Q c h 2 and h, that results 
in the Hubble parameter h being poorly constrained by the 
CMB data alone. This degeneracy is discussed in detail in 
the next Section. 

We note that nearly all of the likelihood is contained 
well within our prior regions, except for the case of tensor 
models with CMB-only data in Fig |^: here there is a region 
allowed by CMB outside our priors with high tensor frac- 
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Figure 3. As Fig. ^| but now considering a wider class of models that possibly include a tensor component. 



tion, h > 1, 



1.3, Q. c h 2 ~ 0.06. These parameters are 



ruled out by many observations apart from 2dFGRS, so the 
truncation is not a concern. 



3.2 Results 

The recovered mean and standard rms error calculated for 
each parameter (except r which is effectively unconstrained) 
are given in Table ^. What is striking is how well specified 
many of the parameters are. 

The general features are as follows: changing from the 
Wang et al. compilation to our compilation slightly shrinks 
the error bars (due to VSA and CBI), but the central values 
are similar except for a slight shift in n s . Allowing tensors 
widens the error bars and causes modest shifts in central 
values (the best fit has a zero tensor fraction, but the fact 
that r must be non-negative explains the shifts). The CMB 
data alone constrains Q,bh 2 and n s well and Q c h 2 quite well, 
but Q,rn and h less well. Adding the 2dFGRS data shrinks 



the errors on fi c /i 2 , h and thus Q, m and fib /fi m by more than 
a factor of 2. 

The most restrictive case is the set of scalar-only mod- 
els. These yield h = 0.665 with only a 7% error, which 
is substantially better than any other method. The mat- 
ter density parameter comes out at Q m = 0.313, with a 
rather larger error of 18%; errors on h and $l m are anticor- 
related so the physical matter density is well determined, 
tt m h 2 = 0.136 ± 7%. We show below in Section ^ that this 
is because the CMB data measure very accurately the com- 
bination Qmh 3 , so that an accurate measurement of Q m re- 
quires h to be known almost exactly. 

Moving from matter content to the fluctuation spec- 
trum, the scalar-only results give a tantalizing hint of red 
tilt, with n s — 0.963 ±0.042. Current data are thus within a 
factor of 2 of the precision necessary to detect plausible de- 
grees of tilt (e.g. n s — 0.95 for \<j> 4 inflation; see Section 8.3 
of Liddle & Lyth 2000). Inflation of course cautions against 
ignoring tensors, but it would be a great step forward to rule 
out an n s — 1 scalar-only model. 
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Table 4. The recovered mean and root mean square (rms) error for each parameter, calculated by marginalizing over the remaining 
parameters. Results are presented for scalar-only and scalar +tensor models, and for CMB data only or CMB & 2dFGRS power spectrum 
data. To reduce round-off error, means and rms errors are quoted to an accuracy such that the rms error has 2 significant figures. We 
also present constraints on some of the possible derived parameter combinations. (Note that due to the marginalization, the maximum- 
likelihood values of 'derived' parameters e.g. f2 m are not simply ratios of the ML values for each 'independent' parameter) 





parameter 


results: scalar only 


results: with tensor component 






CMB 


CMB±2dFGRS 


CMB 


CMB±2dFGRS 




tt b h 2 


n none; _u n nnoo 
U.UzUO ± U.UUzz 


U.UZiU ± U.UUzl 


n noon _i_ n nnQi 
U. Lizzy ± U.UUol 


U.UzZo ± U.UUzo 




fl c h 2 


0.118 ± 0.022 


0.1151 ±0.0091 


0.100 ±0.023 


0.1096 ±0.0092 


using the 
CMB data 
compilation 
of Section |^ 


h 


U.u4 ± U.1U 


U.ODO ± U.U4/ 


U. f O ± U.lo 


n 7nn _L n hkq 
U. f UU ± U.Uoo 


n s 


nnrrn I c\ r\ a a 

0.950 ± 0.044 


0.963 ± 0.042 


1.040 ± 0.084 


1.033 ± 0.066 


m 
r 






n no _L_ o i a 
0.09 ± 0.16 

0.32 ± 0.23 


c\ no _l_ n 1 a 
0.09 ± 0.16 

f\ q o _l_ n oo 
0.32 ± U.ZZ 




0.38 ± 0.18 


0.313 ±0.055 


0.25 ±0.15 


0.275 ±0.050 






0.226 ± 0.069 


0.206 ± 0.023 


0.174 ±0.063 


0.190 ±0.022 






0.139 ± 0.022 


0.1361 ±0.0096 


0.123 ±0.022 


0.1322 ±0.0093 






0.152 ± 0.031 


0.155 ±0.016 


0.193 ±0.048 


0.172 ±0.021 




n b h 2 


0.0209 ± 0.0022 


0.0216 ±0.0021 


0.0233 ±0.0032 


0.0233 ± 0.0025 




Ci c h 2 


0.124 ±0.024 


0.1140 ±0.0088 


0.107 ±0.025 


0.1091 ±0.0089 


using the 
Wang et al. (2002) 
compilation 


h 


0.64 ±0.11 


0.682 ±0.046 


0.74 ±0.14 


0.719 ±0.054 


n s 


0.987 ± 0.047 


1.004 ±0.047 


1.073 ±0.087 


1.079 ± 0.073 


nt 
r 






0.03 ±0.15 
0.25 ±0.21 


0.03 ±0.15 
0.27 ±0.20 






0.41 ± 0.20 


0.296 ±0.051 


0.28 ±0.17 


0.261 ±0.048 






0.240 ± 0.076 


0.200 ±0.021 


0.189 ±0.071 


0.185 ±0.021 




n m h 2 


0.145 ± 0.024 


0.1356 ±0.0092 


0.131 ±0.024 


0.1324 ±0.0088 




Hj, /fi m 


0.149 ± 0.033 


0.160 ±0.016 


0.186 ±0.049 


0.177 ±0.021 



Including the possibility of tensors changes these con- 
clusions only moderately. The errors on h and Q. m hardly 
alter, whereas the error on n s rises to 0.066. The preferred 
model has r = 0, although this is rather poorly constrained. 
Marginalizing over the other parameters, we obtain a 95% 
confidence upper limit of r < 0.7. One way of ruling out the 
upper end of this range may be to note that such tensor- 
dominated models predict a rather low normalization for 
the present-day mass fluctuations, as we now discuss. 



3.3 Normalization 

An advantage of the new CMB data included here is that 
the most recent experiments have a rather small calibration 
uncertainty. It is therefore possible to obtain precise val- 
ues for the overall normalization of the power spectrum. As 
usual, we take this to be specified by as, the rms density 
contrast averaged over spheres of 8 h~ x Mpc radius. For the 
scalar-only grid of models shown in Fig. 2, this yields 

ag = (0.72 ± 0.03 ± 0.02) exp r. (3) 

The first error figure is the 'theory error': the uncertainty in 
as that arises because the conversion between the observed 
Ci and the present P(k) depends on the uncertain values of 
Qm etc. The second error figure represents the uncertainty 
in the normalization of the Ci data (see Fig. 7). The total 
error in as is the sum in quadrature of these two figures. 

This result confirms with greater precision our previous 
conclusions that allowed scalar-only models prefer a rela- 
tively low normalization (Efstathiou et al. 2002; Lahav et 
al. 2002). As discussed by Lahav et al. (2002), a figure of 
as ~ 0.72 is consistent with the relatively wide range of 



estimates from the abundance of rich clusters, but is lower 
than the as — 0.9 for Sl m ~ 0.3 preferred by weak lens- 
ing studies. The obvious way to reconcile these figures is 
via the degenerate dependence of as on r. The lowest plau- 
sible value for this is r = 0.05, corresponding to reioniza- 
tion at z r = 8 for the parameters given here. To achieve 
as ~ 0.9 requires r = 0.22, or reionization at z r = 22, which 
is somewhat higher than conventional estimates (z r < 15; 
see e.g. Loeb & Barkana 2001). Additional evidence in this 
direction comes from the possible first detection of Sunyaev- 
Zeldovich anisotropies at £ > 200 by the CBI (Mason et al. 
2002). This signal is claimed to require as — 1 (Bond et al. 
2002), which would raise z r to almost 30. Further scrutiny 
of these independent estimates for as will be required before 
one can claim evidence for first object formation at extreme 
redshifts, but this is an exciting possibility. 

Finally, we note that these problems are sharpened if 
the CMB power spectrum has a substantial tensor compo- 
nent. As shown by Efstathiou et al. (2002), the model with 
the maximal allowed tensor fraction (r — 0.6) has a nor- 
malization lower by a factor 1.18 than the best scalar-only 
model. This pushes z r to almost 40 for as = 1, which starts 
to become implausibly high, even given the large uncertain- 
ties associated with the modelling of reionization. 



4 THE HORIZON ANGLE DEGENERACY 

In this section we explore the degeneracy observed in 
Figs. |^ & ^ between £l c h 2 and h. This is related (but not 
identical) to the geometrical degeneracy that exists when 
non-flat models are considered, and we now show that it is 
very closely related to the location of the acoustic peaks. 
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Below, we first review the basics of the geometrical degen- 
eracy, secondly note why this is only weakly broken by the 
flatness assumption, and thirdly give a simple heuristic argu- 
ment why this degeneracy approximately follows a contour 
of nearly constant f2 m /i 3 . 



4.1 The geometrical degeneracy 

The 'geometrical degeneracy' in the CMB is well known 
(Zaldarriaga et al. 1997; Bond et al. 1997; Efstathiou & Bond 
1999). If we take a family of models with fixed initial per- 
turbation spectra, fixed physical densities uj m = Q 
Qf,h 2 , and vary both Qa and the curvature Qk to keep a fixed 
value of the angular size distance to last scattering, then the 
resulting CMB power spectra are identical (except for the 
integrated Sachs-Wolfe effect at low multipoles which is hid- 
den in cosmic variance, and second-order effects at high £). 
This degeneracy occurs because the physical densities w m ,u)b 
control the structure of the perturbations in physical Mpc 
at last scattering, while curvature and A (plus u> m ) govern 
the proportionality between length at last scattering and 
observed angle. Note that h is a 'derived' parameter in the 
above set, via h — [w m /(l — Qk — Oa)] 0,5 , so the geometrical 
degeneracy is broken by an external measurement of h. 



4.2 The flat-universe case 

Assuming a flat universe, Ofc = 0, thus also breaks the ge- 
ometrical degeneracy. However, as noted by e.g. Efstathiou 
& Bond (1999), and investigated below, there is a closely 
related degeneracy based on varying two free parameters 
(chosen from Q m , ui m , h, Qa) so as to almost preserve the 
locations of the first few CMB acoustic peaks. This is illus- 
trated in Fig [l|, where the likelihood contours in the (Q m , h) 
plane for CMB-only data form a long and narrow 'banana' 
with its long axis at approximately constant Q m h 3 . The ba- 
nana is surprisingly narrow in the other direction; this means 
that Q m h 3 is determined to about 12% (la) by the CMB 
data. 

This 'banana' is similar in form to the line in Fig. 4 of 
Efstathiou & Bond (1999), though different in detail because 
they used simulations with u> m = 0.25. It is also similar to 
that in Fig. 1 of Knox et al. (2001) as expected. However 
both those previous papers presented the degeneracy in the 
(Q,A,h) plane; although this is just a mirror-image of the 
(Q m ,h) plane, it is less intuitive (e.g. changing Qa alters 
observables that have no explicit A-dependence, via the con- 
straint Q m = 1 — Qa), so the simple Q m h s dependence has 
not been widely recognised. 



4.3 Peak locations and the sound horizon 

The locations £ m of the first few CMB acoustic peaks may 
be conveniently expressed (e.g. Hu et al. 2001, Knox et al. 
2001) as 



£ m — £a(tti - (fim), m = l,2,3 

£ A = 7V/9 S 

_ r s (z„) 



(4) 
(5) 

(6) 



where r s is the sound horizon size at last scattering (redshift 
z»), D A is the angular diameter distance to last scattering, 
therefore 8 S is the 'sound horizon angle' and £ A is the 'acous- 
tic scale'. For any given model, the CMB peak locations 
£ m (m — 1,2,3) are given by numerical computation, and 
then Eq. (jjj) defines the empirical 'phase shift' parameters 
(j> m . Hu et al. (2001) show that the m 's are weakly depen- 
dent on cosmological parameters and <j>i ~ 0.27, <j>2 ~ 0.24, 
4>-i ~ 0.35. Extensive calculations of the (p m 's are given by 
Doran & Lilley (2002). 

Therefore, although 9 S is not directly observable, it is 
very simple to compute and very tightly related to the peak 
locations, hence its use below. Knox et al. (2001) note a 
'coincidence' that 6 S is tightly correlated with the age of 
the universe for flat models with parameters near the 'con- 
cordance' values, and use this to obtain an accurate age 
estimate assuming flatness. 

4.4 A heuristic explanation 

Here we provide a simple heuristic explanation for why 9 S 
and hence the £ m 's are primarily dependent on the param- 
eter combination Q m h iA . 

Of the four 'FRW parameters Q m , u m , h, Qa, only 2 are 
independent for flat models, and we can clearly choose any 
pair except for (Q m , Qa)- The standard choice in CMB anal- 
yses is (u) m , Qa) while for non-CMB work the usual choice is 
(Q m , h). However in the following we take u) m and Q m to be 
the independent parameters (thus h, Qa are derived); this 
looks unnatural but separates more clearly the low-redshift 
effect of Q m from the pre-recombination effect of u m . We 
take ujt as fixed unless otherwise specified (its effect here is 
small) . 

We first note that the present-day horizon size for flat 
models is well approximated by (Vittorio & Silk 1985) 



7'H 



> = o) = fn 



; ' 4 = 6000 Mpc 



(7) 



(The distance to last scattering is ~ 2% smaller than the 
above due to the finite redshift of last scattering). Therefore, 
if we increase Q m while keeping Lu m fixed, the shape and 
relative heights of the CMB peaks are preserved but the 
peaks move slowly rightwards (increasing i) proportional to 
Qm 1 (Equivalently, the Efstathiou-Bond 1Z parameter for flat 
models is well approximated by 1.94 Qm 1 )- 

This slow variation of £a oc Q m X at fixed u! m explains 
why the geometrical degeneracy is only weakly broken by 
the restriction to flat models: a substantial change in Q m at 
fixed uj m moves the peaks only slightly, so a small change in 
u> m can alter the sound horizon length r s (z,) and bring the 
peaks back to their previous angular locations with only a 
small change in relative heights. We now give a simplified 
argument for the dependence of r s on ui m - 

The comoving sound horizon size at last scattering is 
defined by (e.g. Hu & Sugiyama 1995) 



r s (z*) = 



HoQm 2 J,, 



(a + a aq )^ 



da 



(8) 



where vacuum energy is neglected at these high redshifts; 
the expansion factor a = (l + z)~ x and a,, a eq are the values 
at recombination and matter-radiation equality respectively. 
Thus r s depends on u m and Lib in several ways: 
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(i) The expansion rate in the denominator depends on u) m 
via a eq . 

(ii) The sound speed Cs depends on the baryon/photon 
ratio via Cs = c/y/3(l + R), R = 30496 

(iii) The recombination redshift z* depends on both the 
baryon and matter densities U!b,u>m in a fairly complex way. 

Since we are interested mainly in the derivatives of rs 
with Lo m , cut, it turns out that (i) above is the dominant ef- 
fect. The dependence (iii) of z» on uj m ,ujb is slow. Concern- 
ing (ii), for baryon densities u>b — 0.02, Cs declines smoothly 
from cjy/Z at high redshift to 0.80c/\/3 at recombination. 
Therefore to a reasonable approximation we may take a fixed 
'average' c s — 0.90c/\/3 outside the integral in Eq. (^), and 
take a fixed z t , giving the approximation 



r s (z t ) 



0.90 



m(z = 1100) 



(9) 



where r H is the light horizon size; this approximation is 
very accurate for all ui m considered here and Lib — 0.02. 
For other baryon densities, multiplying the rhs of Eq. 
by (wi,/0.02)~ 07 is a refinement. (Around the concordance 
value u>b = 0.02, effects (ii) and (iii) partly cancel, because 
increasing uib lowers the sound speed but also delays recom- 
bination i.e. increases a*). 

From above, the (light) horizon size at recombination is 



r H (z») = 



1/2 



6000 Mpc 



(a + a^) 1 /'- 



■ da 



yl + (a cq /a») — yj a eq /a* 



Dividing by D A ~ 0.98 r H (z = 0) from Eq. 
angle subtended today by the light horizon, 



1.02 



VTT. 



1 + 




a* y a* 

Inserting z, = 1100 and a cq = (23900 a; m ) _1 , we have 
1.02 a-, 01 



(10) 

gives the 
(11) 



vrToi 



1+0.313 



/ 0.147 X 



0.313 



0.147 



(12) 



and 8 S ~ 8 H X 0.9/\/3 from Eq. This remarkably simple 
result captures most of the parameter dependence of CMB 
peak locations within flat ACDM models. Note that the 
square bracket in Eq. (|l^) tends (slowly) to 1 for a cq <C a„ 
i.e. uj m 3> 0.046; thus it is the fact that matter domination 
was not much earlier than recombination which leads to the 
significant dependence of # H on u> m and hence h. 

Differentiating Eq. ( |l^ ) near a fiducial ui m = 0.147 gives 

d In 6 



d In 0„ 



ain6» H 



= -0.1, 



dhu 



-1/2 



= +0.24, 



(13) 



and the same for derivatives of In 6 S from the approximation 
above. In terms of (O m , h) this gives 



ain( 



d In 57 r , 



+0.14, 



ainf 



d In ft, 



+0.48, 



(14) 




Figure 4. Likelihood contours for Q m against h for scalar only 
models, plotted as in Fig. ^. Variables were changed from f!j,/i 2 
and f2 c /i 2 to f2 m and f2{,/f2 m , and a uniform prior was assumed 
for Q{,/f2 m covering the same region as the original grid. The 
extent of the grid is shown by the dotted lines. The dot-dash line 
follows the locus of models through the likelihood maximum with 
constant Q m h 3 i . The solid line is a fit to the likelihood valley 
and shows the locus of models with constant f2 m h 3 (see text for 
details). 



in good agreement with the numerical derivatives of £ A in 
Eq. (A15) of Hu et al. (2001). Note also the sign difference 
between the two d/d\n£l m values above. 

Thus for moderate variations from a 'fiducial' model, 
the CMB peak locations scale approximately as @-m 
n^ ' 14 h~ 0AS , i.e. the condition for constant CMB peak loca- 
tion is well approximated as Q, m h iA = constant. This con- 
dition can also be written cj m f2~ ' 41 = constant, and we see 
that, along such a contour, u> m varies as 1)J, 41 , and hence 
the peak heights are slowly varying and the overall CMB 
power spectrum is also slowly varying. 

There are four approximations used for 9 S above: one 
in Eq. (Q), two (constant c s and z*) in Eq. and finally 
the Q,mJi A line is a first-order (in log) approximation to 
a contour of constant Eq. (|l2|). Checking against numerical 
results, we find that each of these causes up to 1% error in 9 S , 
but they partly cancel: the exact value of S varies by < 0.5% 
along the contour h = 0.7 (O m /0.3)- 1/3,4 between 0.1 < 
fl m < 1. The peak heights shift the numerical degeneracy 
by more than this (see below), so the error is unimportant. 

A line of constant S7 m /i 3,4 is compared to the likelihood 
surface recovered from the CMB data in Fig. ^ In order to 
calculate the required likelihoods, we have made a change 



of variables from ujb & £l c h to fi m & f2j,/fl m , and have 
marginalized over the baryon fraction assuming a uniform 
prior in f2f,/O m covering the limits of the grid used. As ex- 
pected, the degeneracy observed when fitting the CMB data 
alone is close to a contour of constant £ A hence constant # H - 
However, information about the peak heights does alter this 
degeneracy slightly; the relative peak heights are preserved 
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Multipole j2 

Figure 5. The top panel shows three scalar only model CMB angular power spectra with the same apparent horizon angle, compared 
to the data of Table hi Although these models have approximately the same value of Q m /i 34 , they are distinguishable by peak heights. 
Such additional constraints alter the degeneracy observed in Fig. 4 slightly from Q m h 34 to Q m h 3 . Three scalar only models that lie in 
the likelihood ridge with f2 m /i 3,0 are compared with the data in the bottom panel. For all of the models shown, parameters other than 
Q m and h have been adjusted to their maximum-likelihood positions. 



at constant u) m , hence the actual likelihood ridge is a 'com- 
promise' between constant peak location (constant Q m h 3A ) 
and constant relative heights (constant Q. m h 2 )\ the peak lo- 
cations have more weight in this compromise, leading to a 
likelihood ridge along Q, m h 3 ~ const. This is shown by the 
solid line in Fig. ^j. To demonstrate where this alteration is 
coming from, we have plotted three scalar only models in the 
top panel of Fig. |H| These models lie approximately along 
the line of constant Q. m h 3A , with fi m = 0.93,0.36,0.10. Pa- 
rameters other than fi m and h have been adjusted to fit to 
the CMB data. The differing peak heights (especially the 3rd 
peak) between the models are clear (though not large) and 
the data therefore offer an additional constraint that slightly 
alters the observed degeneracy. The bottom panel of Fig. g 
shows three models that lie along the observed degeneracy, 
again with Q m — 0.93,0.36,0.11. The narrow angle of in- 
tersection between contours of constant Q m h 3A and fl m h 2 
(only 10 degrees in the (lnfi m ,ln/i) plane) explains why the 
likelihood banana is long. 

The exponent of h for constant # H varies slowly from 



2.9 to 4.1 as u> m varies from 0.06 to 0.26. Note that Hu 
et al. (2001) quote an exponent of 3.8 for constant l\\ the 
difference from 3.4 is mainly due to the slight dependence 
of 4>\ on ujra which we ignored above. However since that 
paper, improved CMB data has better revealed the 2nd and 
3rd peaks, and the exponent of 3.4 is more appropriate for 
preserving the mean location of the first 3 peaks. Also, note 
that the near-integer exponent of 3.0 for the likelihood ridge 
in Fig. ^ is a coincidence that depends on the observed value 
of Lu m , details of the CMB error bars etc. However, the ar- 
guments above are fairly generic, so we anticipate that any 
CMB dataset covering the first few peaks should (assuming 
flatness) give a likelihood ridge elongated along a contour of 
constant Q m h p , with p fairly close to 3. 

To summarise this section, the CMB peak locations are 
closely related to the angle subtended by the sound horizon 
at recombination, which we showed is a near-constant frac- 
tion of the light horizon angle given in Eq. (|l^). We have 
thus called this the 'horizon angle degeneracy' which has 
more physical content than the alternative 'peak location 
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h 

Figure 6. Likelihood contours for the equation of state of the vac- 
uum energy parameter w against the Hubble constant h. Dashed 
contours are for CMB+2dFGRS data, solid contours include also 
the HST key project constraint. Contours correspond to changes 
in the likelihood from the maximum of 2Aln£ = 2.3,6.0,9.2. 
Formally, this results in a 95% confidence limit of to < —0.52. 



degeneracy'. A contour of constant 8 S is very well approxi- 
mated by a line of constant Q m h 3A , and information on the 
peak heights slightly 'rotates' the measured likelihood ridge 
near to a contour of constant Q m h 3 ' . 



5 CONSTRAINING QUINTESSENCE 

There has been recent interest in a possible extension of 
the standard cosmological model that allows the equation 
of state of the vacuum energy w = p V ac/pva C c to have 
w 7^ — 1 (e.g. Zlatev, Wang & Steinhardt 1999), thereby not 
being a 'cosmological constant', but a dynamically evolv- 
ing component. In this section we extend our analysis to 
constrain w; we assume w does not vary with time since a 
model with time-varying w generally looks very similar to 
a model with suitably-chosen constant w (e.g. Kujat et al. 
2002). The shapes of the CMB and matter power spectra 
are invariant to changes of w (assuming vacuum energy was 
negligible before recombination): the only significant effect 
is to alter the angular diameter distance to last scattering, 
and move the angles at which the acoustic peaks are seen. 
For flat models, a useful approximation to the present day 
horizon size is given by 



r H (z = 0) = -^fi 
-no 



-2w 



(15) 

1 - 3.8w v ; 

(compare with Eq. [7] for w = — 1). As discussed previously, 
the primary constraint from the CMB data is on the angle 
subtended today by the light horizon (given for w = — 1 
models by Eq. ^2|). If w is increased from —1 at fixed 57 m , h, 
the peaks in the CMB angular power spectrum move to 
larger angles. To continue to fit the CMB data, we must 
decrease Q. m h ZA to bring 6 H back to its 'best-fit' value. 



However, the 2dFGRS power spectrum constraint limits 
n m h = 0.20±0.03, so to continue to fit both CMB+2dFGRS 
we must reduce h. 

The CMB and 2dFGRS datasets alone therefore con- 
strain a combination of w and h, but not both separately. 
The dashed lines in Fig. ^ show likelihood contours for w 
against h fitting to both the CMB and 2dFGRS power spec- 
tra showing this effect. Here, we have marginalized over Q m 
assuming a uniform prior 0.0 < Q m < 1.3. An extra con- 
straint on h can be converted into a limit on w: if we in- 
clude the measurement h = 0.72 ± 0.08 from the HST key 
project (Freedman et al. 2001) we obtain the solid likeli- 
hood contours in Fig. []. The combination of these three 
data sets then gives w < —0.52 (95% confidence); the limit 
of the range considered, w — —1.0, provides the smallest 
uncertainty. The 95% confidence limit is comparable to the 
w < —0.55 obtained from the supernova Hubble diagram 
plus flatness (Garnavich et al. 1998). See also Efstathiou 
(1999), who obtained w < —0.6 from a semi-independent 
analysis combining CMB and supernovae (again assuming 
flatness) . 



6 PREDICTING THE CMB POWER 
SPECTRUM 

An interesting aspect of this analysis is that the current 
CMB data are rather inaccurate for 20 < £ < 100, and 
yet the allowed CDM models are strongly constrained. We 
therefore consider how well this model-dependent determi- 
nation of the CMB power spectrum is defined, in order to see 
how easily future data could test the basic CDM+flatness 
paradigm. 

Using our grid of ~ 2 x 10 s models, we have integrated 
the CMB+2dFGRS likelihood over the range of parameters 
presented in Table |^ in order to determine the mean and rms 
CMB power at each i. These data are presented in Table ^ 
at selected I values, and the range of spectra is shown by the 
grey shaded region in the top panel of Fig. |^. A possible ten- 
sor component was included in this analysis, although this 
has a relatively minor effect, increasing the errors slightly (as 
expected when new parameters are introduced), but hardly 
affecting the mean values. The predictions are remarkably 
tight: this is partly because combining the peak-location 
constraint on Q m /i 3 ' 4 with the 2dFGRS constraint on Q m h 
gives a better constraint on Q. m h 2 than the CMB data alone, 
and this helps to constrain the predicted peak heights. 

The bottom panel of Fig. ^ shows the errors on an 
expanded scale, compared with the cosmic variance limit 
and the predicted errors for the MAP experiment (Page 
2002). This comparison shows that, while MAP will beat 
our present knowledge of the CMB angular power spectrum 
for all I < 800, this will be particularly apparent around the 
first peak. As an example of the issues at stake, the scalar- 
only models predict that the location of the first CMB peak 
should be at I — 221.8 ±2.4. Significant deviations observed 
by MAP from such predictions will imply that some aspect 
of this model (or the data used to constrain it) is wrong. 
Conversely, if the observations of MAP are consistent with 
this band, then this will be strong evidence in favour of the 
model. 
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Figure 7. Upper panel: the grey shaded region shows our prediction of the CMB angular power spectrum with 1<t errors (see text); 
points show the data of Table |lj The lower panel shows fractional errors: points are the current data, dashed line is the errors on our 
prediction, and the three solid lines are expected errors for the MAP experiment (Page 2002) for Al = 50 and the 6 month, 2 year and 
4 year data (top-bottom). The dotted line shows the expected cosmic variance, again for A£ = 50, assuming full sky coverage (the MAP 
errors assume 80% coverage). As can be seen, the present CMB and LSS data provide a strong prediction over the full Grange covered 
by MAP. 



7 SUMMARY AND DISCUSSION 

Following recent releases of CMB angular power spectrum 
measurements from VSA and CBI, we have produced a new 
compilation of data that estimates the true power spectrum 
at a number of nodes, assuming that the power spectrum be- 
haves smoothly between the nodes. The best-fit values are 
not convolved with a window function, although they are 
not indep endent. The data and Hessian ma trix are avail- 
able from ittp : //www. roe . ac .uk/~wjp/CMB/. We have used 
these data to constrain a uniform grid of ~ 2 x 10 s fiat 
cosmological models in 7 parameters jointly with 2dFGRS 
large scale structure data. By fully marginalizing over the 
remaining parameters we have obtained constraints on each, 
for the cases of CMB data alone, and CMB+2dFGRS data. 
The primary results of this paper are the resulting param- 
eter constraints, particularly the tight constraints on h and 
the matter density Q m : combining the 2dFGRS power spec- 
trum data of Percival et al. (2001) with the CMB data 



compilation of Section ^| we find h = 0.665 ± 0.047 and 
Q m — 0.313 ± 0.055 (standard rms errors), for scalar-only 
models, or h = 0.700 ± 0.053 and ft m = 0.275 ± 0.050, al- 
lowing a possible tensor component. 

We have also discussed in detail how these parameter 
constraints arise. Constraining Q to t = 1 does not fully break 
the geometrical degeneracy present when considering mod- 
els with varying fitot, and models with CMB power spectra 
that peak at the same angular position remain difficult to 
distinguish using CMB data alone. A simple derivation of 
this degeneracy was presented, and models with constant 
peak locations were shown to closely follow lines of constant 
Q m h 3A . We can note a number of interesting phenomeno- 
logical points from this analysis: 

(i) The narrow CMB f2 m — h likelihood ridge in Fig. ^ 
derives primarily from the peak locations, therefore it is in- 
sensitive to many of the parameters affecting peak heights, 
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Table 5. The predicted mean and rms CMB power calculated 
by integrating the CMB+2dFGRS likelihood over the range of 
parameters presented in Table pi allowing for a possible ten- 
sor component. These data form a testable prediction of the 
CDM+flatncss paradigm. 
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e.g. tensors, n s , r, calibration uncertainties etc. Of course it 
is strongly dependent on the flatness assumption. 

(ii) This simple picture is broken in detail as the cur- 
rent CMB data obviously place additional constraints on the 
peak heights. This changes the degeneracy slightly, leading 
to a likelihood ridge near constant Q m h s . 

(iii) The high power of h 3 means that adding an external 
h constraint is not very powerful in constraining Q m , but an 
external Q m constraint gives strong constraints on h. A 10% 
measurement of Q m (which may be achievable e.g. from evo- 
lution of cluster abundances) would give a 4% measurement 
of h. 

(iv) When combined with the 2dF power spectrum shape 
(which mainly constrains Q m h), the CMB+2dFGRS data 
gives a constraint on fi m /i 2 = 0.1322 ± 0.0093 (including 
tensors) or f2 m /i 2 = 0.1361 ± 0.0096 (scalars only), which is 
considerably tighter from the CMB alone. Subtracting the 
baryons gives Q c h 2 = 0.1096 ± 0.0092 (including tensors) 
or Q, c h 2 — 0.1151 ± 0.0091 (scalars only), accurate results 
that may be valuable in constraining the parameter space 
of particle dark matter models and thus predicting rates for 
direct-detection experiments. 

(v) We can understand the solid contours in Figure ^ sim- 
ply as follows: the CMB constraint can be approximated as 
a 1-dimensional stripe Qmh 3 - = 0.0904 ± 0.0092 (including 
tensors) or f2,„/i 3 ' = 0.0876±0.0085 (scalars only), and the 
2dF constraint as another stripe Q. m h = 0.20 ± 0.03. Mul- 
tiplying two Gaussians with the above parameters gives a 
result that looks quite similar to the fully-marginalized con- 



tours. In fact, modelling the CMB constraint simply using 
the location of the peaks to give Q, m h 3A = 0.081 ± 0.012 
(including tensors) or Q, m h iA — 0.073 ±0.010 (scalars only) 
also produces a similar result, demonstrating that the pri- 
mary constraint of the CMB data in the (Q m , h) plane is on 
the apparent horizon angle. 

In principle, accurate non-CMB measurements of both 
Q m and h can give a robust prediction of the peak locations 
assuming flatness. If the observed peak locations are signifi- 
cantly different, this would give evidence for either non-zero 
curvature, quintessence with m / -1 or some more exotic 
failure of the model. Using the CMB data to constrain the 
horizon angle, and 2dFGRS data to constrain Q m h, there 
remains a degeneracy between w and h. This can be bro- 
ken by an additional constraint on h; using h = 0.72 ± 0.08 
from the HST key project (Freedman et al. 2001), we find 
w < —0.52 at 95% confidence. This result is comparable to 
that found by Efstathiou (1999) who combined the super- 
novae sample of Perlmutter et al. (1999) with CMB data to 
find w < —0.6. 

In Section ^| we considered the constraints that com- 
bining the CMB and 2dFGRS data place on the CMB an- 
gular power spectrum. This was compared with the pre- 
dicted errors from the MAP satellite in order to determine 
where MAP will improve on the present data and provide 
the strongest constraints on the cosmological model. It will 
be fascinating to see whether MAP rejects these predictions, 
thus requiring a more complex cosmological model than the 
simplest flat CDM-dominated universe. 

Finally, we announce the public release of the 2dFGRS 
power spectrum data and associated covariance matrix de- 
termined by Percival et al. (2001). We also provide code 
for the numerical calculation of the convolved power spec- 
trum and a window matrix for the fast calculation of the 
convolved power spectrum at the data values. The data are 



available from either pttp : //www. roe . ac .uk/~wjp/ 


tittp : // www . mso . ami . edu . au/ 2dFGRS 


; as we have 



strated, they are a critical resource for constraining cosmo- 
logical models. 
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